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We present a simple argument to show that the /3-function of the d- 
dimensional KPZ-equation (d > 2) is to all orders in perturbation 
theory given by 



"qq"' Kay Jorg WieseQ 
<3\ 

I—* 
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Neither the dynamical exponent z nor the roughness-exponent £ have 
any correction in any order of perturbation theory. This shows that 
c/2 \ standard perturbation theory cannot attain the strong-coupling regime 

■ and in addition breaks down at d = 4. We also calculate a class of 

correlation-functions exactly. 
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During the last years, there has been an increasing interest in out of equilibrium dynamics. 



Among these, a lot of research was devoted to non-linear growth, and in particular to the 
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Kardar-Parisi- Zhang equation 



GV dh(x,t) _ 2 A ^ 



vV 2 h(x, t) + - (Vh(x, t)y + r](x, t) , (1.1) 



dt y ' ' 2 



77(x, t)r/(a; / , f ) = 2 J D5 d (x - - t') . (1.2) 



Thanks to a fluctuation dissipation theorem and the mapping to exactly solvable models, 



much is known for space-dimension d = 1 |T], 0, |3[]. In contrast, the case of d > 2 can only 
be attacked by approximative methods or field-theoretic perturbative expansions. Using 
the latter, the fixed point structure of the renormalization group flow for d = 2 + e has 
been obtained JI], |], |!|. Two domains can be distinguished: For small effective coupling 

£ = — 3-, (1.3) 



/y 3 



the renormalization group flow goes to in the long- wavelength limit. For large coupling 
the flow is expected to tend to a strong coupling fixed point g = g sc . The crossover takes 
place at g = g co , which turns out to be of order e in an e-expansion and can therefore be 
studied perturbatively. 
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In this article we present a simple argument to resum the perturbation expansion 
and to calculate the renormalization group functions to all orders. This topic has first 
been adressed in §, but is difficult to access there by a non-specialist. The author of 
the present publication was therefore encouraged to find a simple derivation, which sets 
the results on a clear footing and allows to study the limits of the method. Emphasis is 
laid upon a pedagogical presentation, understandable with an elementary background in 
renormalization theory. We will therefore perform all steps of the renormalization program 
by using elementary tools only. This includes a proof of perturbative renormalizability, 
which in standard field theories is a formidable task, see e.g. M and references therein. 

Let us also mention that similar conclusions have independently been obtained by 
H.K. Janssen ||, and in a different context and with completely different methods in 



2 Summation of the KPZ-equation to all orders in perturbation 

theory 



First of all, we want to eliminate the nonlinear term in Eq. (|1 . 1|) . Using the well-known 
Cole Hopf transformation 

W(x,t) :=e^ x ti , (2.1) 
and absorbing a factor of v into t leads to the following equations in terms of W(x, t): 



%-W(x,t) = AW(x,t) + 2L v (x,t)W(x,t) 
at 2u z 



rj(x, t)r)(x', t') = 2uD5 d {x - x')5{t - t') 



(2.2) 
(2.3) 



In interpreting Eqs. (|1 . 1|) and ( [2.2| ) in Ito-discretization, we have explicitly subtracted a 
drift term ~ W(x, t). Thus the expectaton value of W(x,t) will be constant. 

We are now in a position to write down the generating functional for the dynamic 
expectation values 



V 



W 



V [W] V [rj\ e - J K^] + / j(x,t)w(x,t)+j(x,t)w(x,t) 



(2.4) 



with 



J 



W,W,7] 



x,t 



W(x, t) ( W(x, t)-AW(x, t) - -^r](x, t)W(x, t) 



T] 2 (x,t) 



(2.5) 



Expectation values are obtained from Eq. ( |2.4p through variation with respect to j(x,t) 
and j(x,t). The interpretation of the functional ( |2.5|) is simple: The term proportional 
to W(x, t) is just the equation of motion Q2.2p , thus integrating over all purely imaginary 
fields forces the equation of motion to be satisfied. The last term in Eq. ( |2.5| ) is the noise 
distribution. Note that the path integral runs over positive values of W(x, t) only, since 
W(x,t) = e^*'*). 
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The noise-integration can be done. We obtain a simplified action 

.'/ < - ■ ■ ^ 



J 



W, W 



W(x,t) (\V(x,t) - AW{x,t)) - ~ (w(x,t)W(x,t)] 

(2.6) 

where 

9 = —- ( 2 - 7 ) 

As a side remark, let us note that another way to obtain Eq. ( |2.6[ ) is to write down the 
generating functional for the original KPZ-equation (1.1) and then to perform a change 
of coordinates || 

W(x,t) :=e^ h ^ , 

W(x,t) := h(x,t)e-£ h(x ' t) . (2.8) 

This transformation leaves the integration measure invariant. 

Eq. (|2.2j) only makes sense when specifying the initial conditions, i.e. W(x,t) at time 
t = 0. The simplest choice W(x, 0) = leads to d t W(x, 0) = and consequently to 
W(x,t) = 0. We therefore start with 

W(x,0) = l, (2.9) 

which is equivalent to a flat initial condition for h(x,t), namely h(x, 0) = 0. In order to 
eliminate the constant part of Eq. ( |2.9|) from perturbation theory, we set 

W{x,t) = l + w(x,t) . (2.10) 

The response-function of the non-interacting theory ( "free response-function" ) is 

R{x -x',t- t') = (w(x, t)W(x', t')) 

= Q(t - t') [4vr(t - t')Y d/2 e -(— ') 2 /4(t-t') _ (2.11) 

All other free expectation values vanish 

(w(x,t)w(0,0)) o = , 
(W(x,t)W(0,0)) = 0. (2.12) 

Let us now adress the problem of restricting the path-integral to values of W(x, t) > 0. 
Starting with W(x,t) > 0, the time evolution in Eq. (|2.2| ) will keep W(x,t) > for all 
t. This is easily verified for vanishing noise, therefore the free response-function ( 2.11 ) is 
correct. We shall see below that also perturbation theory respects this property. 

Perturbation theory is developed by starting from the functional fl2.6|) . The non-linear 
term is 

~ (W{x, t)W{x, t)f = \-W 2 [x, t) + W 2 {x, t)w{x, t) + i (W(x, t)w(x, t)f , (2.13) 
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and we denote 



1 



1 
2 

1 



2 Jx,t 



W 2 (x,t) = C , 

x,t 

W 2 (x,t)w(x,t) = -C , (2.14) 
(W(x,t)w(x,t)) 2 =^x: ■ 



Since < and - C by its own can not build up divergent diagrams, we neglect them for the 
moment and start by analysing the perturbative expansion of an observable O with ~Xi 
only 

(O) = (Oe*X) . (2.15) 
The basic ingredient is the exponential of the interaction 

e 9 ^ , (2.16) 

from which we have to build vertices in perturbation theory. First of all, there is no 
vacuum-correction, as self-contractions of >C vanish identically due to causality. This 
also holds for the contraction of more than one vertex. With the same argument, we 
conclude, that no diagram with two external legs can be constructed. Therefore, there 
is no divergent contribution to both WW = Ww and W-^Aiy = WAw at any order in 
perturbation-theory, therefore v (hidden in t) has not to be renormalizedQ. The only 
possible diagrams are chains of >C , of the form "."X7X , and so on or higher 

order vertices. The latter are irrelevant in perturbation theory jTJ. 
We therefore write 

1 j — 'V — '\, — '\v 1 j — / \» — '\> — '\, — '\> 1 

+ higher order vertices , (2-17) 

where the time- argument of the vertices grows from left to right. Note that the combi- 
natorial factor of which comes from the expansion of the exponential function at order 
g n has canceled against the n! possibilities to order the vertices in time. In addition, any 
bubble appears with a combinatorial factor of 2, which cancels against factors of 1/2 from 
the vertex, Eq. Q2.13) . So any of the chain diagrams in ( p. 16 ) still contains a factor of 
1/2- 

To proceed further, we first suppress the "higher order vertices" in Eq. (|2.17|) , as the 
only divergencies they may contain are sub-chains as those depicted in Eq. ( p. 17] ) , that 
will be treated here. 

Second, we can switch to Fourier-representation, thus regard the diagrams in Eq. ( |2.17| ) 
as a function of the external momentum p and frequency uo instead of the coordinates x 

2 Note that this is not in contradiction with the non-trivial value for z obtained in dimension d = 1: 
there, as well as for d > 2, the non-trivial fixed point describing the rough phase is in the strong 
coupling regime, i.e. not acceccible by a systematic perturbation expansion; this means that the expansion 
parameter is always large and the expansion uncontrolled. 
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and t, and finally integrate over p and ui instead of x and t. Then, each chain in Eq. 
( [2.17| ) factorizes, i.e. can be written as product of the vertex >C times a power of the 
elementary loop diagram (which is a function of p and uj) 



~yc~yc~yc^x^x' = (o 5 / ~yc ■ (2.18) 



n loops 

Eq. fl2.17| ) is a geometric sum, equivalent to 



1 .'/ 

and one reads off the effective 4-point function 



^ + 9-, — : IX . ( 2 - 19 ) 



^wwWW 



= 9 , • (2-20) 



As we shall show below, the loop integral in Eq. ( |2.20| ) is divergent for any p and u when 
d — > 2. Renormalisation means to absorb this divergence into a reparametrization of the 
coupling constant g: We claim that there is a function a = a(d), such that the 4-point 
function is finite (renormalized) as a function of g R instead of g, when setting 

g = Z g g R p,~ £ (2.21) 

with 



Z g = - , e = d-2 . (2.22) 



1 

1 + ag R 

fi is an arbitrary scale, the so-called renormalization scale. As a function of g R , the 4-point 
function reads 

9RH~ e 



r 
1 



1 + (a-n £ Op^Js 
To complete the proof, we have to calculate the elementary diagram, 



This is 



2 f- fc 



W) * ft + *)»+,(, + „) (|_,) 2 + i ( f _.) 



;2.23) 



(2.24) 



(2.25) 



To perform the integration over i/, the integration path can be closed either in the upper 
or lower half-plane. Closing it in the upper half-plane, we obtain: 

d d k 1 r°° , r d d k 



ds / - — — e 



W d (* + §) +(*-§) +iw J ° J (2*0 

= -J_ r dss -d/2 e - s( i p2+ H 

(8vr) d / 2 7o 



1 /l 2 N^-^A eT 



-p' + iuj) r M - - I . (2.26) 
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The 4-point function in Eq. (p. 23 ) therefore depends on p and u, and more specifically on 



the combination |p 2 + We now chose a subtraction scale, i.e. we demand that T wwWW 
evaluated at /i 2 = \p 2 + iui be 



^wwWW 



9r • (2.27) 



This is achieved by setting 



(8vr) d / 2 V 2 J (8n) d / 

Moreover, since 



r 1 - « = 7Z 2 - - - . (2.28) 



1 /l \d/2-l 

-p 2 + iu) ii~ £ (2.29) 



9r • (2-30) 



e V2 J 

is finite in the limit s — > as long as the combination of |p 2 + icu is finite, it can be read 
off from Eq. ( |2.23| ) that then also ^ ww y/w ls finite. (If useful, either p = or uj = 

may safely be taken.) This completes the proof. Note that this ensures that the model 
is renormalizable to all orders in perturbation-theory, what is normally a formidable task 
to show [J/J. 

The /^-function that we shall calculate now is exact to all orders in perturbation theory. 
As usual, it is defined as the variation of the renormalized coupling constant, keeping the 
bare one fixed 

(3(g R ) = ^ ^ 

From Eq. ( |2.27| ) we see that it gives the dependence of the 4-point function on p and u> 
for fixed bare coupling. Solving 

9 = T— 2.31 

1 + ag R 

for gn, we obtain 

9R = -zr » (2-32) 

fi £ -ag 

and hence 

P(9r) = eg R {l + ag R ) . (2.33) 
Using a from Eq. (|2.28| ), our final result is 




P(jg R ) = {d- 2)g R - -J—J.T 2 - - \ g R , (2.34) 



as stated in the abstract. It shows that standard perturbation theory fails to produce a 
strong coupling fixed point, a result which cannot be overemphasized. This means that 
any treatment of the strong coupling regime has to rely on non-perturbative methods. It 
does of course not rule out the possibility to find an exactly solvable model, non-equivalent 
to KPZ, for which it is possible to expand towards the strong-coupling regime of KPZ. 
Note that also for d — 1 this equation does not possess a fixed point describing the 
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rough phase; the latter is in the strong coupling regime, not accessible by a perturbation 
expansion. 

Let us also note that the /3-function is divergent at d — 4, and therefore our perturba- 
tion expansion breaks down at d — 4. To cure the problem, a lattice regularized version 
of Eq. ( |2.2| ) may be used. However, then the lattice cut-off a will enter into the equations 
and the result is no longer model- independent. This may be interpreted as d = 4 being 
the upper critical dimension of KPZ, or as sign for a simple technical problem. Compare 



also 1 1 



3 Critical exponents and connection to calculations in the standard 

representation of the KPZ-equation 

Quite a lot of work has been done by directly working in the KPZ-picture. A crucial point 
to understand is therefore the relation of the KPZ-picture and the Cole-Hopf transformed 
model used here. The relation is easy on the level of the /9-functions, which are the same 
for both models. In the Cole-Hopf picture it is also easy to see that, as there are no 
corrections to the response function, v and by this means the dynamical exponent z does 
not acquire any perturbative correction, and thus 

z = 2 (3.1) 

to all orders. 

What is not so easy to compare are correlation functions of the height- variable h(x, t) 
in the original KPZ-language with objects of the Cole-Hopf transformed theory. One has 
to compare, see Eq. (|2.1|) 

(h(x, t) h(x', t')) ~ (In W(x, t) In W(x', t')) . (3.2) 

As logarithms are difficult to handle, one can also study expectation values of vertex- 
operators 

with arbitrary k and p. Expectation values of h can be reconstructed by using the identity 



In W = lim - 



1 - w~ s ~ 


— lim 




e~ K - e~ KW 




5^0 


J K 





(3.4) 



o 



Therefore the characteristic functions contain much more information than e.g. the 2-point 
function 

(W(x,t)W(x',t')} . (3.5) 

Using the tilt-invariance of the KPZ-model ||12|| , Eq. ( |3.1| ) also implies that the rough- 
ness exponent £h, defined via (for ( h < 0; for ( h > one would use (^(h(x,t) — h(x',t)) 2 
instead) 

(h(x, t) h(x', t)) ~ \x - x'\ 2(h (3.6) 
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of the /i-field at a non trivial fixed point is 

Ch = . (3.7) 



To study correlation-functions more directly, we use the first term from Eq. Q2.14Q to 
calculate 

r x t 
("•(•'•• ') «•(•'''•/')) .</ / y,rC ' 

Jy,r x ,t 



gf R(x-y,t-T)R(x'-y,t'-r) . {3.1 
•>y,r 



Note that this is the only diagram which contributes, since more complicated diagrams 
involving loops like <d j-XZ. > nave to be taken at zero momentum and frequency and thus 
vanish according to Eq. ( |2.26 ). 



Further, when t and t' are small, the expectation value in Eq. Q3.8|) is small, too. This 
is physical, since the surface has not much grown yet. In the other limit of large times, 
oo and keeping t — t' fixed, the r.h.s. of Eq. ( |3.8|) converges towards 

gC(x-x',t-t') , (3.9) 

where C(x,t) is the standard dynamic correlation function which in Fourier space reads 

C ^) = TT—2- ( 3 - 10 ) 
For equal times, this relation reads 

(w(x, t) w(x', t)) ~g\x- x'\ 2 ~ d = g R Z g \x - x'\ 2 ~ d , (3.11) 
leading to a renormalization for w of the form 

w R = Zg*W (3.12) 
and to a roughness exponent ( w for the field w in 

(w(x, t) w(x', t)) ~ \x - x'\ Kw (3.13) 

with 

(w = Co + $(w i 

Co = — , (3-14) 

Using Eqs. ( |2.21| ) and ( |2.22| ), we can express (3(g) in terms of 5( w (gji) as 

P(9R) = 9R[e-28C w ] ■ (3.15) 
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At a non-trivial fixed point gn = g* > 0, i.e. at the roughning transistion, this induces 

SCM = | , (3.16) 

and therefore 

Cw = • (3.17) 

It is tempting to indentify ( w with (h. A priori, this is surprising, since h and w are very 
different observables, related by 

w (x,t) =e$ h(x ^ - 1 . (3.18) 

The putative identity ( w = (h will therefore only hold, if w is, within correlation functions, 
well approximated by the term linear in h on the r.h.s. of Eq. (|3.18| ). It will certainly 
break down when Ch > 0, i.e. in the strong-coupling regime. It is still possible to relate 
correlation-functions for h and w via integral-transforms as in Eq. (|3.4|) , as long as the 
expectation value in Eq. ( |3.3j ) is dominated by contractions with only, leading to purely 
Gaussian correlations. However, in the strong-coupling regime, also non-linear terms, of 
which the first is (w 2 (x,t)w(x',t')), will contribute to Eq. (|3~3| ), making a systematic 
analysis easier said than done. This is another difficulty of the perturbation expansion 
beyond the roughning transition. 

On a more technical level, it is worth realizing that the above relations can be used 
to simplify the perturbation expansion in the original KPZ-language. 



4 Conclusions 

In this article, we have presented a simple method to resum the perturbation expansion 
of the KPZ-equation and to calculate the renormalization group /3-function to all orders 
in perturbation theory, including a proof of perturbative renormalizability. The main 
conclusion is that there is no anomalous contribution to the dynamical exponent z in 
the weak-coupling regime and at the roughening transition. We also have given some 
indications of why standard perturbation theory fails in describing the strong-coupling 
fixed regime. 
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